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ABSTRACT 


Earthquake data from Central Calatfornia cVe reer oes 
January 1, 1969, to December 31, 19/715 and £romiicereun 
Wevdada over the period Auplist 31, 1954, tomMecenber = 1 on, 
have been analyzed to investigate the hypothesis that the 
Farth eves may be the triggering mechanism for yeiewoceut penec 
©f earthquakes. One method used compared the tidal components 
and the component rates-of-change at the time and location of 
Meclidlwearengquakes to those for random events Unitvonnl, 
Beneradeedmover “awe Same time period. A second methem con- 
Paned the drrterences in the Earth tides computed Gdeross tic 
Ereaath Of a fault zone at the time of aGtual ecarthe@makes fo 
mee Ger terences at the times of random events unitorml, 
Belewaeed over the same time period. Statistical tests were 
performed to see if the distributions from the actual events 
were the same as the corresponding distributions from the 


randomly generated events. 
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I. INTRODUCTION 


Every point on the surface of the Eawth 1s subjeceuce 
taree types of forces: the force of eravity ducetomene 
attraction of the mass of the Earth, the centrifugal force 
ame to the rotation of the Earth, and a force due to the 
attraction of extra-terrestrial bodies, mainly the Sun and 
tne Moon. The resultant forces on the surface of the Earth 
mececenerally referred to as the Earth tide forecs. 

Because of the stresses that these Earth tides Grearcuum 
mae crust of the Earth, it has been long felt that theyemay 
be a significant factor in the triggering of earthquakes. 
several previous studies have been undertaken in which 
meremptsS were made to correlate the gravitational compoment 
of Che Earth tides with the occurrence of earthquakes. In 
one study, Knopoff [Ref. 1] concluded that there was no 
correlation between this tidal component and the occurrence 
of the 3800 earthquakes in Southern California that he inves- 
tigated. Shlien [Ref. 2] came to a similar conclusion in his 
study of several widely separated regions of the world. 

A team of Russian researchers [Ref. Si eontend that enece 
and other similar studies have suffered from two serious 
drawbacks. First, they considered only the gravitational 
component of the Earth tides. More important may be the 
tangential components which depend very strongly on the 
Orientation of a fault in a region.  In@oeneral etre extreme. 
in the tangential components do not occur at the times of the 


extremes of the gravitational component. Second, they tried 
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to cover too large an aréa. Mitepehtseonsidered arco ace 
18 square degrees, and the regions im Shilien ses tid) ea eo 
between 50 and 450 square degrees. Since tides are time 
dependent and periodic in nature, considerations of large 
memrons have averaging effects on their variation » possi, 


masking dependencies present. 





Il. OBMSCTIVE 


The objective of this study was three-fold Gm aed 
region of investigation and a specific time period: (1) 
ECompare the tidal forces from a series of actual Carthenakes 
to these forces from a series of random events generated 
uniformly over the same time period, (2) compare the rates- 
of-change of the tidal forces from the actual events to those 
from the randomly generated events, and (3) Compare the 
differences in each of “the tidal”f6rces over the™breadth of 
the region from the actual events to those from the randomly 
Memweraved events. These comparisons were used €o tesemunc 
Ingeernesis that the distribution of a foreey rate ,eeuneaiere r 
ence from the actual events was the same as the corresponding 
distributions from the randomly generated events. Rejection 
@f the hypothesis wotld indicate a difference in the distri- 
butions being compared. Thus, an evaluation can be made as 
to whether a significant dependency exists between a force, 
rate, or difference and the occurrence of earthquakes for the 


particular region being studied. 





III. ORGANIZATION 


For the purpose of thiS analysis, the ti@altG@me oon wee 
resolved into five quantities. The orientation and con- 
Struction of these quantities were=as follons- 9) tjeetia a 
tional component was formed along the radius of the Earth at 
mene C€arthquake epicenter, the normal Component was tommpea 
perpendicular to both the gravitational component and the 
Meerike of a fault, the axial component was formed perpendic- 
ular to the gravitational component and parallel to the 
Strike of a fault, the tangent magnitude was the resultant 
formed by the combination of the magnitudes of the normal 
and axial components, and the total magnitude was the 
mesultant formed by the combination of the magnitudes of the 
gravitational, normal, and axial components. By observing 
Mmewethese five quantities changed over a smal) changewm 
mime, the rate-of-change of each quantity also was appne a. 
mated. By calculating these tidal quantities on both sides 
of a fault zone, the differences In Ghese Gquantitresmas © 
were calculated. 

It should be noted here that the@tonegome det in rere momen 
are mormal and axial components of thewtraal forces usm 
peccludes their being exactly tangemtmtoutace surtace ormene 
Barth Owing to the fact that the @banriiienered peuces 
Sphere. osimee this angular devilatien fren tangency was so 
Small and was nearly constant over the small range in lati- 
tude of the regions investigated, no correction was mMaqdevin 


Calculating component values. This correction should eames 





out during the comparisons of component frequency 


distributions and, therefore, have little influence on the 
meswilts of the comparisons. 

Invsiellecting regvons forMinvestigatron, Zones@e. hips 
seusmac sactivity were desined so that rela@mwely small regions 
mould be selected and stall yield a sufficient numbervor 
earthquake occurrences to support a detailed analysis. The 
mefcions chosen were small sections of the San Andreas taulet 
aneCentral California and the Dixie Valley fault in Western 
Nevada, both of which are located in the Western United 
States (Figure 1). These two regions were chosen because of 
maear differing fault structure. The San Andreas fault @isoa 
Serike-Slip fault where the hand mass to the southwest side 
of the fault moves northward in relation to the land mass on 
mpageiOomtheast side of the fault (Figure 2). In the partic- 
ter region of the San Andreas fault investigated, the strike 
@iethe fault 1s approximately 40 degrees west of north as 
measured ren longitude 121 W. The Dixte Valley fauieeree 
graben fault where the material in the center of the fault 
Zenec moves downward in relation to the sides of the fault 
(Figure 3). The strike of this fault is approximately 
mecallel to longitude 118 W. Co-ording@tes of these regio. 
ame given in Table 1. 

i= CiOestsicman LIMe per Lodmaven which earthquakes were to 
be analyzed, a large enough time period was needed so that 
bias due to the periodic fluctuation or suncmear chetaces 


(monthly, annually, etc.) would be damped or averaged out. 





This was needed so that differences found in the Comparison 
Hf these quantities would not be attributable Vomtnewe ime 
mature of the tides. For the specified portion Of eee om 
Andreas fault, 1732 earthquakes of Richter magnitude 0.0 or 
Preater occurred in the three year period from Januar, 1. 
ieee to Vecember 51, 19715 For the por fron ores sD 
Valley fault, the available data prior to 1962 was for shocks 
®iL magnitude 4.0 or greater, and from 1962 on, for shocks of 
Magnitude 0.0 or greater. So as not to unduly bias the 
ometri bution of data with the increased frequency of recorded 
small shocks after 1962, only the time period through 1961 
was considered. During the seven year period from August 31, 
1954 (the start of detailed recording in this area), to 
December 31, 1961, 100 earthquakes of magnitude 4.0 or greater 


occurred. 
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IV. EARIEA Dore Cus orton 


The actual tidal data were generated through the Wsiemoea 
Eomputer routine (unpublished) by Professor R. He shuddee 
Naval Postgraduate School. Much work has been accomplished 
Ser the past few years on the subject of theoretical Earen 
tides. Melchior [Ref. 4] has one of the better discussions 
@metnie subject. He goes into great detail concerning the 
eomputation of the theoretical Earth tides using an equi- 
potential surface argument. His procedure determines the 
ii@eadtions of the Sun and the Moon in relation to the Earth ae 
meearcrwcular time and establishes the potential surftaceceior 
M@m@estides. tIhen by specifying the location of a point on ythe 
Earth in relation to these surfaces, he determines the tides 
at that location. Shudde's program takes a different approach. 
imealso determines the locations of the Sun and the Moon usage 
eiemprocedures im reference 5S, Dut then at ae particulan spon 
@omene Larth, the tides are calculated directly using tic 
standard gravitational force laws. For a more detailed 
description of the Earth tide calculations, see the discussion 
in Perendix N. Pt will suffice Nerve suey, that all Shudde'‘s 
program requires for inputs are the time, date, and geographic 
co-ordinates of an earthquake occurrence. For this analysis, 
the times and locations for the actual occurrences were taken 
from references 6-8 and from data supplied by the National 
Genter for Earthquake Research. The dates @ion ene Tandem 


occurrences were generated uniformly over the time periods 
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specified in Section II. Them@ii@ention for these aan 
events was taken to be the center of the reanonmeeeune 


mimes tigated. 
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Ve ANALYSIS 


The objective of this analysis was to test them goocn-. 
ae ao 
maere A is the distribution of data from the actualecycmeae 
and R is the distribution of the corresponding data from the 
randomly generated events. Rejection of this hypothesis for 
Maeparticular force, rate, or difference would indicate a 
difference in the two distributions being compared and, hence, 
acependency between that quantity and the occurrence of 
earthquakes. 
ie primary statistical test used to test this hypoemeanr- 
was the Kolmogorov-Smirnov (K-S), Two-Sample test. This is a 
non-parametric test independent of distributional assumptions 
meemt the data other than the distribution be continuous. | Et 
is an exact test and operates by the comparison with the 
actual data rather than with grouped data. A second test 
Meta Was the Chi-Square Goodness-of-Fit test which operaces 
by the comparison of grouped data. 
The K-S Mavis ic "nossesses the advantage that it is an 
exact method, whereas the Chi-Square method requires a 
fairly large sample to justify the approximations that are 
meeded sto applyewit'' [Hoel, 1971]. 
The Chi-Square test was used because this is the common 
statistical test used in many studies. The K-S test was used 
because it is a more powerful test and was more appropriate 


to the data. 


A. KOLMOGOROV-SMIRNOV, TWO-SAMPLE TEST 


This test was implemented through the use of the IBM 
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computer subroutines KOLM2Z [Refs 10], with a moditited scan 
routime [Ref. lb]. Usang these routines, for eachwof ithe 

five tidal, rate, and difference quantities, thewdistyibwmaon 
of data from the actual earthquakes in a region was compared 
tomehes distribution for that quantity calcullatedstpome€ he 
nmandemly generated events. The two sets of data for each 
distribution were each sorted into increasing order and formed 
into sample distribution functions. The maximum absolute 
S@mainate difference between the two distributions was deter- 
mined. A sample size adjustment factor was calculated by 
fmecane the square root of the ratio of the sample size 

fjmacauiet to the sample size sum. The K-5, two-sample staereiic. 
Z, Was then calculated by taking the product of the maximum 
difference and the adjustment factor. That is, 


Z = max |A(x)- S Goff a: a 


where A(x) is the distribution cane pee Beite actual eveneee 








S(x) is the distribution function from the random events, Ny 
1s the number of events in A(x), and N. is the number of 
@yents in S(x). The use of the adjustment factor 1s necessany 
in the K-S test. If not applied to the statistic as shown 
above, the eee of the factor would need to be applied to 
the critical values listed in Table 2 to keep the test statis- 
tic consistent with the critical values. In this case, appli- 
cation of the adjustment factor to the statistic simplified 
the use of the computer routines and enabled the critical 


values to remain constant for a given level of significance. 


ey 





Be CHI-SQUARE GOODNESS -OF [3 leew 

For each tidal quantity, the maximum range covering Doth 
the tidal data from the actual earthquakes and the data from 
mire randomly génerated eventsS were divided inte 20 equim 
macervals and the data was grouped into these intervals. 
mapacent intervals were collected together unt) these peer 
meeauency an the collected amtervalsecxeccdc amen ce moo elce 


W965}. The Chi-Square statistic was then calculated as 
m! 2 
= Ss” (0; -E:)} 
“¢ a 
b=] 
mere O. is the frequency in interval} 1, trom the actwal 
feuts, E; is the expected frequency Gm interval, 1, angen 


meethe numbereef collected intervals. Selected critical 


marues for the Chi-Square test are given in Table 3. 
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mop ted from Beyer, p. 294. 


Critical Values ot ene 
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Table 3 
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VI. RESvEis 


In the analysis of the tidal data, seme bilas Geen 
distributions of the data was expected due to the physical 
misctance between the origin of the actual earthquakes vanderme 
moometric center of the region at which the random evene- 
were generated. Although this bias was felt to be minimal 
mae to the small size of the regaons and the llarse number ver 
Seientcs, additional checks were performed to test the size oF 
this bias. Random events were generated at the north, south, 
Meet, and east extremes of each region and the five tidal 
S@m@e@erive rate quantities were calculated at™each Pocatworm for 
aieme occurrence. The ten quantities at the north and se@men 
extremes were then compared using the K-S test, as were the 
memeauantities from the west and east extremes. Tne results 
are shown in Tables 4 and 5. The conclusion was that bias in 
the data due to the physical separation of the events was so 
slight as to have an insignificant influence on the final 
analysis. It mus t be stressed heremthat the sicenificanecmoen 
the following results pertains only to the particular region 
and pHe nature and orientation of the particular fault being 


analyzed. 


A. SAN ANDREAS 

The K-S, two-sample statistics for each of the five 
mraal, tive rate, and five difference Gquantatics ace 17 em 
the San Andreas region are shown in Table 6 and the Chi- 


Sauare statistics for these quantities are shown in” lables 


Za 





Referring to Tables 2 and 3 for a 0.05 leveliceetico ie ae 
the normal tidal and gravitational ratesquantitiesmi en 
S@atastically significant umder both tests. None of the 
tidal differences was statistically significant at this level 


in either test. 


fe DIXIE VALLEY 

Me K-S, two-sample statistics for each of the tivo ome 
mimemnate, and five difference quantities acting in the wire 
Peiiey region are shown in Table 8 and the Chi-Square statrs-— 
mies tor those quantities are shown in Table 9.  Referrinoeio 
iees 2 and 3 for a 0.05 level of significance, the norma 
rate quantity was statistically Significant under both team, 
Heme Of the tidal difference quantities was significant at 
mies level in either test. 

miter much of this analysis had been completed, a 
Semalarity between the methods of calculating the rate and 
feererence quantities was observed. The rate quantities were 
calculated by computing the tidal quantities at a location, 
allowing a small unit of time to elapse, and then computing 
eS eagle quantities again. The change in the tidal quantity 
divided by the change in time was used as its rate of change. 
Beeause of the rotation of the Earth, changes in time corres= 
pond to changes in location at a fixed latitude. A lapse of 
four minutes is equivalent to a distance change of one degree 
in longitude. For small time lapses, the distance considered 
moved is negligible. For larger time lapses, the calculations 


can be considered as being at two separated locations. The 


bef 





difference quantities were calculated by compucing sige an 
eously the tidal quantities at two locations separated by one 
jiegree in longitude. Since both procédures involved 
mebeulation of components at sepdrated loeati0ns ssereme 
Should have been a correlation between the test results on 
the rate quantities and the test tesults on the difference 
G@uantities. Comparison of the San Andreas results in Table 

6 showed some similarity but the interpretation was clouded 
in that the tangential components in this region were rotated 
mearough an angle of 40 degrees in order to align them with 
om@eestrike of the fault. For the region in which the 
Surzential components had not been rotated, comparison Of the 
temte Valley results in Table 8 showed a marked agreement in 
meemorder of the significance of the results. Therefore, 
ma@esc interpretations of the results could have been made 
from a single analysis of the rate quantities without doing 


mee difference calculations. 
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K-S .~Statiseic 
Component N-S W-E 


Gravitational Tide Re 44 
Normal Tide ig oeo VAs 
Axial Tide “oe 28 
Tangent Magnitude soe 25 
Total Magnitude oe nya. 
Gravitational Tide Rate 42 5 oe 
Normal Tide Rate B72 J 2.0) 
Axial Tide Rate 225 25 
Tangent Magnitude Rate 34 ao 
Total Magnitude Rate 015 ~42 


Results of Bias Check in Tidal Gomearicens 
N-S and W-E Extremes of the San Andreas Area 
Using the K-S  Seacmeene 


Table 4 
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K-5 Gta stic 

Component Neo Wee 
Gravitational Tide oe mad 
Normal Tide . 28 ako 
Axial Tide - 48 . 20 
Tangent Magnitude vole , 16 
Total Magnitude OO AS 
Gravitational Tide Rate 46 sale 
Normal Tide Rate oo ale 
Axial Tide Rate 26 Pees) 
Tangent Magnitude Rate . 30 20 
Total Magnitude Rate . on . 16 


Results of Bias Check in Tidal Comparisons 
N=-S and*WeE Extremes of themii ies) wise aed 
Usins the K-S Staere one 


Table 5 
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Component Statietic 


Gravitational Tide 1.1411 
Normal Tide iodo 
Peatal Tide eee 7 © 
Tangent Magnitude 1.3344 
Total Magnitude 0.8658 
Gravitational Tide Rate e515 
Bewnal Tide Rate OAkZ0 
Peetai lide Rate 0.8469 
Tangent Magnitude Rate | 0.9541 
Total Magnitude Rate Wea 1S 
Geavitatienal Tide Difference Lol 24 
Homemal Iide Difference a 2 SZ 
foamed Tide Difference 0; Sueiey 
Tangent Magnitude Difference lev 2 59 
Total Magnitude Difference i dos 


Results of Actual Tidal Comparisonce mii one 
San Andreas Area Using the K-S Statistic 


Table 6 
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Component 


GSGavitational Tide 
Normal Tide 
Axial Tide 
Tangent Magnitude 


Total Magnitude 


Geavitational Tide Rate 


Normal Tide Rate 
Feral Tide Rate 
Tangent Magnitude Rate 


Total Magnitude Rate 


Seats te 


114. 


40. 
OO: 


36 


Zoe 


34. 


a2 
84 


rg 


ay 


53 


ZA 


-o8 
48 
Gel 
.64 


Degrees of Preedom 


io 
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17 
18 
18 


ie 
18 
Ly 
18 
Jey 


Results of Actual Tidal C@iiparvsoms iu sence 
San Andreas Area Using the Chi-Square Statistic 


Table 


‘eI 
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Component EEE ULSIE TC 
Gravitational Tide ous 
Normal Tide B62 78 
Axial Tide noo 
Tangent Magnitude le 
Total Magnitude moe Se 
Seavitational Tide Rate ~O9 55 
Normal Tide Rate ~4119 
Axial Tide Rate oe Un 
Tangent Magnitude Rate ~4673 
Total Magnitude Rate Ares) 
Mmevwitational Tide Difference ~5743 
Mermal lide Difference 2054 
feral! Tide Difference .6040 
Tangent Magnitude Difference A 
Total Magnitude Difference -Lo0s 


Results of Actual Tidal Comparisons in the 
Dixie Valley Area Using the K-S Statistic 


Table 8 
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Component 


Gravitational Tide 
Normal Tide 
Axial Tide 
Tangent Magnitude 


Total Magnitude 


Gravitational Tide Rate 
Normal Tide Rate 

Axial Tide Rate 

Tangent Magnitude Rate 


Total Magnitude Rate 


Stat iettac 


lig 

6. 
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Loe 


NO 


Lois 


ioe 


Le 
90 
30 
04 
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34 


47 


sale 
sch 


72 


Degreesmot Freedom 


13 
14 
14 
15 


Wa! 


14 
14 
Ws 
EZ 
MS 


Results of Actual Tidal Comparicen- sim scie 
Dixie Valley Area Using the Chi-Square Statistic 


Table 
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VIL. CONCEUSECNS 


Results of the analysis indicate that, even at high lévels 
of significance, there exist components of the Earth tides 
amd their rates-of-change that are Sipnifieantly relaceamnes 
the occurrence of earthquakes and that these components vary 
with the nature and orientation of the fault being investi- 
gated. 

The significance of the normal tidal component on the 
frerike-slip San Andreas fault gives credence to the theory 
jm@at the occurrence of earthquakes in this fault zone may 
mesult from the lessening of the frictional force acting 
Memween the two.vertical faces of the fault. Since the 
fjrrecional force acting yereneerr two surfaces is a function of 
the normal force acting on those surfaces, a lessening of the 
menial force could decrease the frictional force enougiero 
femrow the accumulated stress in the tault to be released: 

ihe relationship between the rate-of-change of a) tiraae 
memec and the difference in a force over the breadth of a 
fault has been discussed previously. Although ie dintewentee 
quantities were considered significantedueaeue os | leve low 
Melfi cance, reference to Tables 2, 6 and 8 reveal that these 
G@m@antities did differ considerably and could have been consid- 
ered Significant at lower levels of sipniticamec” | [he sion, 
meance of the gravitational rate im the regiomeor ones oa 
Andreas fault indicates a differential in the wventical fore. 
acting on opposite sides of the fault.” “timsmsieni: cance 


difficult to interpret in view of the strike-slip nature of 


30 





the fault. Because of the preximity of this automa 
Pacific Ocean, this signifseance may mimdicate ties re ene men 
ecean tidal. loading effects on the ocelrrence cfteameraaa 
methasregion. The significance of the normal ytaceuaue- 
meeion of the Dixie Valley fault indicates a daitfterenera iam 
the east-west force acting on opposite sides of this fault. 
wfias interpretation is consistent with the nature of the 

Dixie Valley fault in that the downward slippage of the center 
peethe fault zone occurs with» the separation of the sidee 


aeche fault. 


om 





APPENDIX A 


A. OBSERVATIONS ON ASSOCIATED LEGENDRE POLYNOMIALS. 


The generating function for the assocraged Wea waa. 


th 
polynomial of the m Side ie 


m 
2 | TE Zz w oOo 0 ra 
= ( ) a mae ad yu Yo. a (x) 
Mn! (i-2zxyy ) Se ja. 


ther useful a aniee 


Py (x) = (1=x) 2 (x) 





3 0OBe Sore, Ther (x) = 2. Py, C)] 


Letting m=l1 in (1) and oe ae terms, 


3 

(\-2xy +7?) ae 2 (ee y* oye ie '(x) 
Prom (2), 

(1 —x?) "Di (y) = | 

(1- x) * P(e) = 3x 

(i = x*) . ee) = (Gm ~1) 

(1-x 2 P(x) = 26x 334) 
From (3), | 

me ee) =| 

im CH) a 


ae (x) = =z (3x*-1) 
P, (x) = +(sx?-3x) 
a (x) - | (a5 x"- 30x +3) 
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(3) 


(4) 


(s) 


(c) 





B. EARTH TIDE CALCULATIONS. 


—— 


S.i0 






fe- Earth radius vector to epicenter 


ee = Scalar magnitude of R 


R,,- = Distance vector from epicenter to B; 

i. — ©calar magnitude of Ieee 

ioe, Dictince VeGctor rom Origin fo Bs 

Bo ~ 2calar magnitude of ae 

ee = Extra-terrestrial body, i=l (moon), 2 (Sun) 
F;, = Force vector between the origin and B, 

= 6 = Force vector between the origin and B; 


%, = Angle between R and R;, 


= 
I 


mass of Bs 


From Newton's laws of gravitation, 


3 k-mi- Rie Ks Mme ° Rs 
pie 7 5 — = 
ree V.o° 


The total tidal force vector, T, is the sum over all bodies 











a T r : k —10 1e 
min 10 Fie i ee as 
LLG ie 
ok 
Since Kem co > 
Ries Sio-R, i = 3 Rio -(Rio-R a ( ) 
oe = ir; pune : 2 





x4 2 ; Lo 
Watt Mie = Tio +a ZaeKio cos (%; ), the expression can 


be written as 
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(3) 





a) |) a(S )eoeor GE) 


Letting x;=cos(4.), and a — ; a comparison of (4) saameh 
Ie 
(8) siete 


(a) (\- la "3" PC) 


Be i. -| 1 P (cos (%: ») 
Foo 


sin (K¢) 





a Q sa A coslhe C4 
=> &) St Py (cos(u:)) ) 


d coslai) 


Ssemstituting (9) into (7) yields, 


_ 4 P (costa: y) 


. Q-| 
Kemi Cc 
: r.? ~ (Rie -R) 2 Kio Hicos On =z ) 
tO 
Letting t, be the component ae Ty in thesdixectionmo.: Niue 
tre vector in the plane normal to R, and letting CP. Dem tne 


angle between N and Ti, wesults in the general expression, 


tyaNl,, and after rearranging terms, 





i = cos( Pi) 5 my oe 5 a (cos(ai )) (10) 


Expanding terms in (10) and substituting values from (5), 


‘.: 7 a cos((;) (c =) + RB (easta) + 








d cos («:) 
Q ) d 72 ( 
rr cos (Ki) ) + ++ 
(F: d cos&:) ) 
Ke mea Q 
=- OES oe(01)| 3 cos(a.)+ 2(--\(sco#@.) - en ( \) 
RR: 
Noting that cos@,)= = =. a 
Saas ae 


Pagieis tdkenmtompe the north(z) dircctrony (1h) egives t, 
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ae 
Tabhecos (ahs — == , and if NaS "taken to bev ine eo Fee 


: 10, eae 
mirection, (11) gives ty with cos(f)= =—: If the 
10 
€o-ordinate system 1S oriented so that Reis sin themx-dviee- 


RT, 
mien, then t x=t, 2 





a 
oS 
= a ig fa. —— L)- Ca: Vio? -cos(& )- a’) +e (Ze) j 
Teeetan Faleos a) | | 


_ kKemi-a@ a | 
=a ales i. costae > - fe om Te hileoss) 
(12) 


momemaing a few terms of (12), collecting coefrieclents io. 





... a score: >= +3(F-\(scos® (Xi) - 3ces(«:))-+ | 


with cos (K;,)= ae For tangential components in other than 
10 


E.) where j=0,1, and substituting values from (6), 





mae y and z directions, N can be any unit vector in the 


meone normal to the vector R. 


Ss 





10. 


ce. 
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